Dynamics of Cold Atoms Crossing a One- Way Barrier 



Jeremy J. Thorn, Elizabeth A. Schoene, Tao Li, and Daniel A. Steck 
Oregon Center for Optics and Department of Physics, 
1274 University of Oregon, Eugene, Oregon 97403-1274, USA 

We implemented an optical one-way potential barrier that allows ultracold *^Rb atoms to transmit 
through when incident on one side of the barrier but reflect from the other. This asymmetric barrier 
is a realization of Maxwell's demon, which can be employed to produce phase-space compression and 
has implications for cooling atoms and molecules not amenable to standard laser-cooling techniques. 
The barrier comprises two focused Gaussian laser beams that intersect the focus of a far-off-resonant 
single-beam optical dipole trap that holds the atoms. The main barrier beam presents a state- 
dependent potential to incident atoms, while the repumping barrier beam optically pumps atoms 
to a trapped state. We investigated the robustness of the barrier asymmetry to changes in the 
barrier-beam separation, the initial atomic potential energy, the intensity of the second beam, and 
the detuning of the first beam. We performed simulations of the atomic dynamics in the presence 
of the barrier, showing that the initial three-dimensional momentum distribution plays a significant 
role, and that light-assisted collisions are likely the dominant loss mechanism. We also carefully 
examined the relationship to Maxwell's demon and explicitly accounted for the apparent decrease 
in entropy for our particular system. 
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I. INTRODUCTION 

The transport dynamics induced by asymmetries of 
systems is an important paradigm in physics, dating 
back to Feynman's thought-experiment of a mechanical 
ratchet at finite temperature [1]. The point of Feyn- 
man's analysis was to demonstrate that in thermal equi- 
librium, the ratchet cannot be used to "rectify" thermal 
fluctuations to do useful work, because Brownian mo- 
tion of the ratchet mechanism itself causes it to intermit- 
tently fail in such a way as to prevent steady-state mo- 
tion against a constant force. However, in nonequilibrium 
situations — such as a simple temperature difference — the 
asymmetry of a system can induce steady-state trans- 
port, and of course do work. Recently, the study of more 
general "ratchets," where thermal fluctuations or time- 
dependent forces couple with periodic asymmetric po- 
tentials to produce steady-state directed motion [2], has 
become of broad interest. Ratchet systems in the form 
of molecular and Brownian motors [3, 4, 5] are of par- 
ticular interest in understanding how controlled motion 
is effected in spite of thermal fluctuations, especially at 
the nanoscale, where such fluctuations tend to be large. 
Ratchets have also recently been studied in the context of 
the laser-cooling of atoms in asymmetric optical-lattice 
potentials [6]. 

A related problem of asymmetric transport occurs for 
particles in the presence of asymmetric barriers. For ex- 
ample, asymmetric diffusion occurs in transport across 
membranes. A recent experiment demonstrated that 
asymmetric transport can be caused by an asymmetry 
in the shape of the membrane pores, provided that a suf- 
flciently wide range of particle sizes is present [7]. In 
this case, larger particles can clog the pores (but only 
from one side), preventing smaller particles from diffus- 
ing through the membrane in one direction. In atom 



optics, it is possible to realize similar "one-way barriers" 
or "atom diodes," as proposed independently by Raizen 
et al. [8] and Ruschhaupt and Muga [9] in slightly differ- 
ent contexts (with a number of subsequent refinements 
in the general design [10, 11, 12]). These optical one-way 
barriers for cold atoms are asymmetric optical potentials 
in the sense that atoms experience a different potential 
depending on the direction of incidence: atoms incident 
from one side reflect from the barrier, while atoms in- 
cident from the other side transmit through it. Unlike 
the above membrane, however, this one-way action is en- 
tirely a single-particle effect, relying on optical potentials 
that change depending on the internal state of the atom. 
Note that in a closely related proposal, a combination 
of coherent and incoherent processes conspires to allow 
transport of particles from one reservoir to another, but 
not in the reverse direction [13]. 

A common theme among all the systems mentioned 
above is that dissipation is required to produce asym- 
metric transport. Ratchets, for example, typically in- 
volve overdamped motion and thus heavy dissipation. 
However, an interesting question is, how little dissipa- 
tion is required to produce a significant asymmetry in 
the transport? The most straightforward answer for the 
optical one-way barrier here is that one photon must be 
spontaneously scattered per particle, since the action of 
the one-way barrier depends on an irreversible change 
in the internal state of the atom. (Presumably, though, 
one can construct more complicated schemes that require 
even less dissipation; we will show below that our one-way 
barrier dissipates more entropy than is required by the 
second law of thermodynamics.) This question is partic- 
ularly relevant in the atom-optical case, since one of the 
main motivations in realizing one-way barriers for atoms 
is in developing new laser-cooling methods for atoms 
[8, 14, 15, 16, 17]. The idea here is that the one-way bar- 
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rier can force the accumulation of atoms into a volume 
much smaller than the original container. The phase- 
space density of the atomic vapor increases if the effect 
of the volume compression outweighs any heating effects 
of the barrier. (Cooling can then be straightforwardly ef- 
fected via adiabatic expansion; the most important step 
in terms of controlling an atomic vapor is reducing the 
phase-space volume.) Phase-space compression with one- 
way barriers has recently been demonstrated with rubid- 
ium atoms [18, 19]. The best increase in phase-space 
density achieved so far is a factor of 350 above the initial 
conditions of a vapor in a magnetic trap [20]. Though 
standard laser-cooHng techniques such as Doppler cool- 
ing are now well established [21], they generally rely on 
the existence of a cycling optical transition: atoms falling 
into "dark states" decouple from the cooling lasers. Since 
the irreversible action of a one-way barrier relies in princi- 
ple only on the scattering of a single photon, laser-cooling 
techniques based on one-way barriers may help circum- 
vent problems with dark states and enable laser cooling 
of molecules [22] and the many atoms for which standard 
laser-cooling methods are ineffective. 

Another, more fundamental, motivation for study- 
ing the one-way barrier comes from its connection to 
Maxwell's demon [15, 23, 24, 25]. In this thought- 
experiment, the demon manipulates a trapdoor in a wall 
dividing a container of gas. The demon could use the 
trapdoor as a one-way barrier to reduce the space oc- 
cupied by the atoms without performing any work, in 
an apparent violation of the second law of thermody- 
namics. (Maxwell's original demon used the trapdoor to 
separate hot and cold atoms to achieve a temperature 
gradient, but the reduction of entropy is the essence of 
the thought-experiment.) One can think about the reso- 
lution to this "paradox" in a number of ways [24, 25] , but 
the key issue is that the demon must perform measure- 
ments on the system, gaining the information needed to 
perform feedback via the trapdoor. The entropy decrease 
of the atoms is balanced by an increase in the entropy 
of the demon's memory due to the accumulated informa- 
tion. This cannot continue indefinitely in a cyclic process 
unless the demon's memory is reset, or "erased." The 
erasure requires transferring entropy to the environment, 
in accordance with the second law. Our experiment ef- 
fectively consists of only a single thermodynamic cycle, 
and so a cyclic erasure is not strictly necessary. The 
spontaneous scattering of photons from the barrier laser 
light in our experiment both acts as an effective position 
measurement on each atom and carries away sufficient 
entropy to compensate for entropy lost in the reduction 
in phase-space volume. 

This work builds on our previous study of an all-optical 
realization of a one-way barrier for rubidium atoms [19]. 
Here we perform a detailed study of the dynamics of cold 
atoms crossing the one-way barrier. We study how var- 
ious aspects of the barrier configuration contribute to 
its operation, and we discuss the key elements of opti- 
mizing the robustness and performance of the barrier. 
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FIG. 1: (Color online) Schematic diagram of the optical setup, 
showing the dipole trap, barrier beams, and imaging system. 
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FIG. 2: Relevant atomic levels of Rb on the D2 cooling 
transition (not to scale), shown with the off-resonant coupling 
of the main barrier beam. Atoms with F=l see this field 
as red-detuned, while F=2 atoms see a blue detuning. The 
repumping barrier beam is resonant with the ^=1 F'—2 
MOT repumping transition. 



We also compare our data to simulations that treat the 
center-of-mass motion of the atoms classically, and ob- 
tain good agreement. Finally, we examine in some detail 
the entropic aspects of phase-space compression with the 
one-way barrier. 



II. PROCEDURES 

The optical one-way barrier is the principal feature of 
our experimental setup (Fig. 1), comprising two focused 
Gaussian laser beams. The one-way barrier scheme ex- 
ploits the hyperfine ground-state structure of ^^Rb to cre- 
ate an asymmetry that allows atoms to transmit through 
the barrier when traveling in one direction, but reflects 
them when traveling in the other. The main barrier beam 
is tuned between the F=l F' and F^2 F' ^'^Rb 
hyperfine transitions, as shown in Fig. 2. Thus the de- 
tuning A :— uj^loq of the barrier-beam frequency lo with 
respect to the atomic resonance ujq has opposite signs for 
atoms in the two hyperfine ground levels. This results in 
the main barrier beam presenting an attractive potential 
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FIG. 3: (Color online) Optical potentials seen by atoms along 
the axis of the dipole trap in our experiment. The overall trap- 
ping potential is from the main dipole-trap beam, while the 
narrow state-dependent feature results from the main barrier 
beam. The main barrier beam produces a repulsive potential 
for atoms in the reflecting {F=2) ground state, and an at- 
tractive potential for atoms in the transmitting (-^=1) ground 
state. 



to atoms in the F=l ground state and a repulsive po- 
tential to atoms in the F=2 ground state, because the 
optical dipole potential is inversely proportional to the 
detuning A. These two potentials are plotted in Fig. 3. 
To create the asymmetry, we introduce a second beam, 
the repumping barrier beam, that is resonant with the 
F=l — > F'=2 repumping transition, displaced from the 
main barrier beam as shown in Fig. 1. If all atoms start 
in the transmitting (F=l) ground state, they can pass 
through the main barrier beam until they pass through 
the repumping barrier beam, at which point they will be 
pumped to the reflecting {F—2) ground state. The atoms 
will then see the main barrier beam as a potential bar- 
rier and will remain trapped on that side of the one-way 
barrier. We choose to have the repumping barrier beam 
on the right-hand side of the main barrier beam, which 
makes that side the reflecting side of the barrier, and the 
left-hand side the transmitting side. 

We initially cool and trap the ^^Rb atoms in a standard 
six-beam magneto-optic trap (MOT) [21], loaded from a 
cold atomic beam produced by a pyramid MOT [26] . Af- 
ter a secondary polarization-gradient cooling stage with 
reduced intensity and increased detuning of the trapping 
beam, we have about 2 x 10*^ atoms at about 30 /iK in 
an ultra-high vacuum of ^ lO"^*' torr. The position of 
the MOT can be shifted several millimeters in each co- 
ordinate using magnetic bias fields. 

After cooling we load the atoms into a far-detuned op- 
tical dipole trap produced by a 1090(5) nm Yb:fiber laser. 
The laser emits a coUimated, multiple-longitudinal-mode, 
unpolarized, nearly Gaussian beam with a 1-9(1) mm 
1/e^ beam radius. We operate this laser so that the total 
power inside the vacuum chamber is 9.3(5) W. We focus 
the beam with a single 200 mm focal length plano-convex 
lens, producing a 30.9(5) /im waist (1 /e^ intensity radius) 
and a 2.8 mm Rayleigh length. For *^Rb atoms in either 
hyperfine ground state, this beam yields a nearly con- 
servative potential well with a maximum potential depth 
of X 0.9 mK. This dipole trap has longitudinal and 
transverse harmonic frequencies of 24 Hz and 3.0 kHz, 



respectively (near the trap center), a 1/e lifetime of 20 s, 
and a maximum scattering rate of only 3 s~^. Atoms are 
typically loaded in the anharmonic region of the trap, 
so that the atomic motion dephases and has a different 
period (50 ms) than the harmonic frequency suggests (an- 
gular momentum adds to this effect; see Appendix C). 

The two one-way barrier beams are nearly parallel 
asymmetric Gaussian beams with a variable separation. 
Their foci nearly coincide with the focus of the dipole- 
trap beam, intersecting it at about 12(3)° from the per- 
pendicular to the beam axis (as in Fig. 1). The main 
[repumping] barrier beam has a waist of 11.5(5) fiia 
[13(2) ^m] along the dipole-trap axis and 80(7) ^m 
[60(7) /im] perpendicular to the dipole-trap axis. We 
control the power of the repumping barrier beam with 
an acousto-optic modulator (AOM). 

We measure the separation of the two beams on a beam 
profiler that has a resolution of 5.6 /im/pixel. We take 
pictures of each beam profile at many locations along the 
beam axis near the focus, thus making it possible to ac- 
curately locate the axial position of the focus despite the 
fact the beam waists are on the same order of magnitude 
as the camera resolution. Using this setup, we determine 
the separation of the beams (on the order of ten microns) 
with an error on the order of a micron. 

We image the atoms by illuminating them with a 45 
pulse of light resonant with the F=2 F'=3 MOT trap- 
ping transition. This absorption-imaging beam is nearly 
perpendicular to the dipole-trap beam, and is detected 
by a charge-coupled-device (CCD) camera, which images 
the shadow left by the atoms that scatter light out of 
the beam. We image atoms in the F=l ground state 
(which do not scatter the MOT trapping light) by turn- 
ing the MOT repumping beams on slightly before and 
throughout the imaging pulse, ensuring that all atoms 
are in the F=2 ground state. The short image pulse is 
ideal for getting accurate spatial imaging, since the atoms 
cannot move much during that time. To reduce system- 
atic errors due to interference fringes in the images, we 
subtract background offsets (as computed from the edge 
regions of the images) on a per-column basis and then 
integrate each column to form distributions (such as the 
ones shown in Fig. 4). The spatial resolution is 24.4 /im 
as set by the CCD pixel spacing, but the distributions 
are smoothed slightly for visual clarity. The total num- 
ber of trapped atoms drifted slowly over time, so we often 
rescaled atom distributions in order to aid comparison of 
different data sets. 

Our "canonical" data are taken as follows. We load the 
dipole trap with the MOT centered on the dipole-beam 
axis 0.95(5) mm away from the dipole-trap focus for 5 ms, 
trapping about 3 x 10'* atoms (measured by imaging the 
resonance fluorescence from the MOT light on a CCD 
camera) at ^100 /iK, with a peak one-dimensional atom 
density of around 4 x lO'' atoms/mm (on the order of 
10^ atoms/mm"^). Longer load times trap more atoms, 
but the atoms spread throughout the trap during load- 
ing. We used a relatively short loading time to keep the 
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atoms localized. To allow atoms from the MOT that 
are not trapped in the dipole trap to fall away, we ex- 
tinguish the MOT beams and wait approximately 20 ms 
(about half an oscillation period) after loading atoms into 
the dipole beam but before turning the barrier beams 
on. During this half-period delay, atoms spread in the 
dipole trap, reducing their apparent mean displacement 
from the dipole-trap focus. For example, the center of 
the atomic distribution for the normal 0.95(5) mm mean 
starting offset of the atoms is only 0.58(8) mm away from 
the dipole-trap focus after the half-period delay. Hence- 
forth, we will quote the offset of the center of the atomic 
distribution after this delay. 

During the half-period delay after loading atoms into 
the dipole trap, we optically pump the loaded atoms to 
the F—1 ground state by extinguishing the MOT re- 
pumping beams and then leaving the MOT trapping 
beams (red-detuned by 70 MHz) on for another 7 ms. 
We verified that this pumps virtually all the atoms to 
the F^l ground state by using absorption imaging with- 
out the MOT repumping beams. 

After the half-period delay, we turn on the barrier 
beams (we designate this time as f = 0), allow the system 
to evolve for some time ranging from a few milliseconds 
to as long as a few seconds, and then image the atoms. 
We load atoms to either side of the barrier beams, and la- 
bel data as either "transmitting side" or "reflecting side," 
based on which side of the barrier beams the atoms were 
on at i = (when the barrier beams are turned on). Be- 
cause we wait a half-period between loading the atoms 
and i = 0, atoms are actually loaded on the opposite 
side of the barrier beams with respect to their "starting" 
location. 

The main barrier beam was normally 34(1) /im to the 
left (as seen by the camera) of the repumping barrier 
beam, making the left-hand (right-hand) side of the bar- 
rier beam the transmitting (reflecting) side. We nor- 
mally ran the main barrier beam with 40(4) fiW of power 
(inside the vacuum chamber), resonant with the *^Rb 
F—3 F'=3, 4 crossover dip in the rubidium saturated- 
absorption spectrum, which is 1.05(5) GHz blue of the 
^''Rb MOT trapping transition. We normally ran the 
repumping barrier beam at 0.36(4) /iW (inside the vac- 
uum chamber), and always stabilized it to the ^^Rb 
F=l F'=2 MOT repumping transition. 

We use two different conventions when reporting laser 
detunings. When we wish to give the absolute frequency 
of a laser, we report the detuning as the difference be- 
tween the laser frequency and the F=2 — > F'=?> *^Rb 
MOT transition. However, for computing optical poten- 
tials and scattering rates, we need to account for the 
presence of four separate excited states, and the hyperfine 
splitting between these states is not negligible. For these 
purposes, we compute and report an effective detuning. 
The dipole potential for each two-level transition is pro- 
portional to the squared dipole matrix element for the 
transition, and inversely proportional to the detuning. In 
order to keep a two-level-atom approximation, we average 



over the different excited states, using the correct dipole 
matrix elements and energy shifts, to obtain an effective 
two- level detuning. For example, in our canonical setup, 
our main barrier beam is tuned 1.05(5) GHz to the blue 
of the F=2 F'=3 *^^Rb MOT transition. When we 
compute an optical potential using all four excited states 
in the D2 manifold to obtain an effective two-level-atom 
detuning that gets the same answer, we find 1.12 GHz for 
atoms in the F=2 ground state, and —5.41 GHz for atoms 
in the F=\ ground state. For computing scattering rates, 
we need to average squared dipole matrix elements di- 
vided by the square of the detuning. For the detunings 
we used, effective two- level detunings for both optical 
potentials and scattering rates differed by less than the 
accuracies quoted here. 

One common variation from our canonical setup was 
the "F=2 setup." For these data, during the half-period 
while we waited for MOT atoms to fall away, we chose 
to optically pump the atoms into the F=2 ground state 
by flashing the MOT repumping beam without the MOT 
trapping beam. Unless otherwise specified, we optically 
pumped atoms into the ^=1 ground state. 

Another variation we used was the "filled trap." For 
these experiments we loaded the atoms to the left-hand 
side of the barrier beams for 110 ms and waited 200 ms 
before turning the barrier beams on, loading about 9 x 
10^ atoms at ~ 100 /iK. These times are several times 
the longitudinal trap period, and so resulted in a fairly 
even distribution of atoms throughout the dipole trap. 
We used the same method described above to ensure all 
atoms were pumped into the _F=1 ground state for these 
experiments. 

Unless otherwise noted, the main barrier beam was 
linearly polarized perpendicular to the dipole-trap beam 
axis, and the repumping barrier beam was linearly po- 
larized parallel to the dipole-trap beam axis. We found 
that beam polarization was not a critical factor in the 
barrier operation. 



III. MAIN RESULTS 

The main results for the one-way barrier are presented 
in Fig. 4, which shows the atomic evolution in response to 
the one-way barrier. To elucidate the effects of the one- 
way barrier. Fig. 4(a) shows the evolution of the atoms in 
the dipole trap in the absence of the barrier. As expected, 
the atoms oscillate back and forth about the trap center 
with some breakup of the atomic cloud evident at later 
times due to the anharmonicity of the trapping potential 
and the spread in angular momentum of atoms orbiting 
the dipole-trap axis. 

Figure 4(b) shows the dynamics of the atoms in the 
presence of the one-way barrier, which is located at the 
origin. The repumping barrier beam is positioned to the 
right of the main barrier beam in the columns so that the 
atoms, initially released from the left-hand side of the 
barrier, will interact first with the main barrier beam. 
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FIG. 5: (Color online) Populations on the left- and right- 
hand sides of the one-way barrier for a "filled-trap" experi- 
ment where atoms were initially distributed throughout the 
dipole trap. Error bars indicate statistical error from 38 rep- 
etitions. 
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FIG. 6: (Color online) Schematic of a four-level atom that 
would allow much larger barrier-beam detunings for reduced 
scattering. Here, the main barrier beam is detuned from a 
different transition than the one with which the repumping 
barrier beam is resonant. If the main-barrier-beam transition 
has a substantially different frequency than the repumping- 
barrier-beam transition, then the main barrier beam may be 
detuned by an amount greater than the ground-state splitting, 
and still create an effective barrier for atoms in the reflecting 
ground state while having negligible effect on atoms in the 
transmitting state. 



The repumping barrier beam defines the reflecting side 
of the barrier, which makes the left side the transmitting 
side. The atoms, which are initially prepared in the F=l 
transmitting state, pass through the main barrier beam 
to the right-hand side of the barrier. There the atoms in- 
teract with the repumping barrier beam, which optically 
pumps them to the F=2 reflecting state. Upon their re- 
turn, the atoms are effectively trapped on the right-hand 
side of the barrier as they repeatedly reflect from the 
barrier. 

When the atoms start on the reflecting side of the bar- 
rier, as in Fig. 4(d), they also remain trapped to the right 
of the barrier. This holds true regardless of the initial 
state of the atoms; atoms starting in the F=2 reflecting 
state are expected to reflect from the main barrier beam, 
while atoms starting in the F~l transmitting state will 
first encounter the repumping barrier beam, which will 
optically pump them to the F=2 reflecting state before 
they reach the main barrier beam. The results presented 
in Fig. 4(d) show the evolution for atoms initially in the 
_F=1 transmitting state, and we observe similarly good 
reflections for atoms starting in the F—2 reflecting state. 

Figure 4(c) reveals what happens when atoms start 
on the left-hand (transmitting) side of the barrier in 
the F=2 (reflecting) state. Initially the atoms reflect 
from the barrier as expected, but then the atoms grad- 
ually pass through the barrier and become trapped on 
the right-hand side. This phenomenon is a direct result 
of choosing the main-barrier-beam frequency to be more 
nearly resonant with the F—2 F' transition than the 
_F=1 — > F' transition. Though the main barrier beam 
reflects many atoms during the initial encounter, this 
interaction will also change the state of many of these 
atoms from the F=2 reflecting state to the F=l trans- 
mitting state. When these atoms encounter the main 
barrier beam the second time, they will transmit through, 
while atoms that did not change state on the first refiec- 



tion will eventually change their state upon subsequent 
refiections. This results in a net accumulation of atoms 
on the right-hand side of the barrier despite the atoms 
being initially in the "wrong" {F—2 refiecting) state. The 
one-way barrier functions almost equally well for atoms 
on either side of the barrier and for atoms in either state. 

As a further test of the utility of the one-way bar- 
rier, we activated the barrier after a "filled-trap" load of 
the dipole trap, which resulted in a symmetric and rel- 
atively uniform initial distribution of atoms, all in the 
F=l (transmitting) ground state. This experiment also 
had a lower barrier-beam power of 18(2) /iW, chosen to 
reduce heating without much harm to the function of 
the barrier. In this scheme atoms on the left-hand side 
of the barrier gradually transmit to the right-hand side, 
while atoms on the right-hand side remain there. The re- 
sults for this test are shown in Fig. 4(e), where many of 
the atoms that start to the left of the barrier eventually 
transmit to the right-hand side. Figure 5 is generated 
from the same data, and shows the populations on each 
side of the barrier as a function of time. Here it is easier 
to see that atoms really are being transferred from the 
left-hand side of the barrier to the right-hand side, in the 
manner of Maxwell's demon. 



A. Scattering 

The main technical challenge associated with imple- 
menting the one-way barrier lies in the limited detunings 
afforded by the 6.8 GHz ground-state hyperfine splitting 
of *^Rb. Detuning the main barrier beam halfway in 
between the F—1 and F=2 ground states seems an ob- 
vious choice; however, this is still near enough to the 
F=l F' resonance to optically pump a portion of 
the atoms (initially in the F—\ transmitting state) to 
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the F=2 reflecting state while they are traversing the 
barrier beam. This results in a sudden increase in po- 
tential energy, greatly increasing the total energy of the 
atom. The heating associated with this transition is on 
the order of the barrier height, which greatly decreases 
the effectiveness of the barrier. This heating mechanism 
informed our choice of main-barrier-beam detuning, in 
which the main barrier beam is more nearly resonant 
with the F=2 F' transition, with a 1.05(5) GHz de- 
tuning from the F—2 F'=3 MOT trapping transition. 
With this detuning, atoms in the F=l ground state scat- 
ter few photons, so heating during transmission is mini- 
mized. Scattering during reflection is more probable be- 
cause the main barrier beam is closer to the resonances 
available to atoms in the F=2 ground state. However, 
when atoms change state during reflection, their poten- 
tial energy is decreased, which cools them. Those atoms 
are then free to transmit, and are still cool enough for 
the barrier to trap them when they return to the reflect- 
ing side. Data collected to probe the effects of different 
main-barrier-beam detunings are presented and analyzed 
in Sec. IV, where we see that a more symmetric detuning 
does not function as well as this asymmetric detuning. 

For our canonical setup, which uses the detuning de- 
scribed above, we expect ~ 0.7 scattering events on a 
single transmission and ~ 8 scattering events on a single 
reflection, ignoring state changes and any scattering from 
the repumping barrier beam (see Appendix B). We also 
ran simulations (see Sec. V for more details) that account 
for state changes, heating, and scattering from the re- 
pumping barrier beam. These simulations found slightly 
higher values of ~ 3 scattering events for transmission 
and ~ 10 scattering events for a single reflection. In par- 
ticular, note that even though the main barrier beam is 
detuned by more than half of the hyperflne ground-state 
splitting for transmitting atoms, scattering is not negli- 
gible. 

Losses in general in this setup seem to be mainly due 
to light-assisted collisions (discussed in Sec. V) facilitated 
by the barrier beams during reflection and transmission. 
Light-assisted collisions include fine-structure-changing 
collisions [27] , hyperfine-changing collisions [28] , photoas- 
sociation [29], and radiative escape [30, 31]. We believe 
that, for our barrier beam, the dominant loss mechanism 
is radiative escape [30, 31], although all can be modeled 
as density-dependent loss mechanisms. In radiative es- 
cape, an atom is promoted to an excited state by one of 
the barrier beams. While excited, the atom is much more 
polarizable and interacts more strongly with neighboring 
atoms. Before decaying, the atom accelerates toward a 
neighboring atom. Once the atom has decayed back to 
the ground state, the atoms cease interacting, but keep 
their kinetic energy. This kinetic energy is often enough 
to eject both atoms from the trap. 

Initial trap lifetimes range from 300 to 500 ms on the 
right-hand (reflecting) side of the barrier depending on 
the temperature of the atoms. The fltted exponential 
lifetimes are markedly larger (approaching 700 to 900 ms) 



in the late-time tails of the population-decay curves. 



The scattering rates are not too problematic in our 
setup, as the lifetimes indicate. However, the lifetimes 
could be improved considerably if we could reduce the 
scattering rate by one or two orders of magnitude, so 
that the average number of scattering events on trans- 
mission and reflection was well below unity. As shown 
in Appendix B, the scattering rate depends mostly on 
the barrier height, atomic speed, and the ratio of the 
linewidth to the detuning. Once the barrier is high 
enough to easily trap atoms, scattering during reflection 
is only weakly dependent on barrier height. Scaling down 
the atomic velocities limits the utility of the barrier as 
a cooling mechanism, so increasing the detuning is per- 
haps the best variable to use in reducing scattering. A 
large increase in detuning could be achieved if the main 
barrier beam used a transition with a substantially dif- 
ferent frequency than any available to the reflecting state 
(see Fig. 6). This would allow the main barrier beam to 
be far-detuned from the reflecting-to-excited state tran- 
sition, and still several orders of magnitude further de- 
tuned from the transmitting-to-intermediate state tran- 
sition. Atoms in the transmitting state would thus be 
largely unaffected by the main barrier beam because of 
the very large detuning, but could still be pumped to the 
reflecting ground state by the repumping beam. If neces- 
sary, a third laser beam detuned from the transmitting- 
to-intermediate state could be used to cancel any shifts 
of the transmitting state by the main barrier beam. This 
would circumvent the ground-state-splitting limitation 
and lead to perhaps orders of magnitude less scattering. 



As a more speciflc example, we could use ^^Sr, with the 
true ground state (5s^ ^So) as the transmitting state and 
the metastable state {5p ^P2) as the reflecting state. For 
the intermediate and excited states shown in Fig. 6, we 
could use 5p ^Pi and 5d "^Ds, respectively. With these 
states, the main-barrier-beam transition is almost 40 nm 
red of the repumping-barrier-beam frequency and has a 
linewidth comparable to the ^''Rb D2 transitions. This 
would allow the main barrier beam to be detuned much 
further than the ^^Rb ground-state splitting and still be 
many nanometers detuned from the repumping-barrier- 
beam transition. This could reduce scattering during 
reflection by multiple orders of magnitude and reduce 
scattering during transmission to negligible levels. The 
repumping decays presented here are slow (millisecond 
time scales) and involve other states that may decay to 
cither the reflecting state or back to the transmitting 
state, but all routes back to the transmitting state are 
short compared with the reflecting state lifetime. The 
long repumping time could be corrected for by allowing 
the repumping barrier beam to flU the entire trapping 
side of the barrier, allowing for more interaction time. 
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IV. ROBUSTNESS 

To investigate the robustness of the one-way barrier, 
we varied four key parameters in the experiment: 

1. the separation between the main barrier beam and 
the repumping barrier beam, from 8 to 74 /im; 

2. the initial displacement of the atoms (this repre- 
sents an effective change in barrier height), from 
to 0.9 mm on either side of the barrier; 

3. the power of the repumping barrier beam, from 
0.002 to 8.7 AiW; and 

4. the detuning of the main barrier beam from the 
F=2 ground state of ^"^Rb, from -1-0.75 to -1-4 GHz 
(blue-detuned) , measured with respect to the *'^Rb 
MOT trapping transition. 

Initially the separation between the beams was set 
to 34 fim, which ensured that the tails of the Gauss- 
ian beam profiles along the dipole-trap axis overlapped 
shghtly. We beheved the overlap was necessary to pre- 
vent the atoms trapped on the right-hand side of the 
barrier in the F—2 reflecting state from being pumped 
back to the F—1 transmitting state by the main barrier 
beam (with an effective two-level A w 27r x 1.12 GHz), 
such that they would pass back through the barrier to 
the transmitting side. The presence of a small amount of 
repumping-barrier-beam light in the same region would 
rapidly pump any atoms that made it to the F—1 trans- 
mitting state back to the reflecting F=2 state. The over- 
lap was designed to guarantee that atoms did not slowly 
leak back across the barrier. 

Surprisingly, our studies of the effect of barrier-beam 
separation indicate that while overlapping the beams is 
important for optimizing barrier performance, it is not 
critical to barrier operation. Figure 7 shows populations 
on either side of the barrier after 100 ms of evolution 
as functions of beam separation. The transmitting-side 
data show the most dramatic effect for separations below 
about 30 /im. When the beams are too close together, 
the tail of the repumping barrier beam can pump atoms 
on either side of the barrier to the reflecting state, caus- 
ing the barrier to reflect atoms from both sides. This 
can be seen for the small separations (on the order of 
two beam waists and below) on the transmitting-side 
data, where the majority of atoms stayed on the left- 
hand (supposedly transmitting) side of the barrier. We 
also see a higher loss rate for small separations, which 
we suspect may be due to atoms changing state in the 
middle of the main barrier beam as opposed to the tail, 
as well as changing state more frequently. When many 
state changes occur in the middle of the barrier beam, the 
large changes in potential energy result in a large amount 
of heating, which causes much more loss. Large separa- 
tions show much less effect on the transmitting-side data, 
although the F=2 data seem to get worse and then im- 
prove as the separation increases past 50 /im. For atoms 



beginning to the right of the barrier, the beam separation 
does not severely affect the barrier's reflectivity, though 
the functionality of the barrier declines as the separation 
is increased beyond about 50 /im. As long as the re- 
pumping barrier beam is still present to the right of the 
main barrier beam, most atoms will pass back through 
the repumping barrier beam after reflection, correcting 
any state changes that may have occurred. As the sep- 
aration between the beams increases, more atoms may 
fail to make it back to the repumping barrier beam af- 
ter multiple reflections, explaining the slow reduction in 
reflectivity beyond 50 /fm. This observation has implica- 
tions for one-way barrier cooling applications that require 
a signiflcant number of reflections from the barrier. 

We investigated how the loading position of the atoms 
in the dipole trap affected the dynamics of the one-way 
barrier as a substitute for changing the barrier height. 
The initial position affects the ratio of average atomic 
kinetic energy to barrier height, and so makes a decent 
substitute for varying the barrier height. However, it is 
important to note that changing the kinetic energy with- 
out adjusting the barrier height is not quite the same as 
changing the barrier height without altering the kinetic 
energy; the two changes have different effects on the ex- 
pected number of scattering events while the atoms tra- 
verse the barrier. Appendix B contains a more detailed 
analysis of this assertion. 

In Fig. 8, we plot the atomic populations in the one- 
way barrier at 100 ms as functions of initial loading po- 
sition, as determined by the location of the center of the 
MOT trapping fields. We see no obvious effect of the 
loading position on the efficiency of the one-way barrier. 
Figure 9 shows some sample time-series for various initial 
loading positions. Essentially, we see that atoms starting 
on the transmitting side are transferred to the reflecting 
side, and atoms starting on the reflecting side are kept 
there, independently of initial loading position. This in- 
dicates that the efficiency of our one-way barrier is fairly 
insensitive to the velocity of the incident atoms. 

While the intensity of the main barrier beam controls 
the height of the barrier, the intensity of the repumping 
barrier beam also plays a crucial role in the transmission 
and reflection of atoms from the barrier. We varied the 
intensity of the repumping barrier beam by over three or- 
ders of magnitude to observe its effect on the functioning 
of the barrier. The results for atoms that started on the 
left-hand and right-hand sides of the barrier are presented 
in Fig. 10, where the populations of atoms on either side 
of the barrier after 100 ms are measured for different 
repumping-barrier-beam intensities. Even over the full 
range of intensities we studied, the barrier continued to 
act asymmetrically in both experiment and simulations. 
This can be seen in Fig. 11, which has time-series plots 
of the highest and lowest repumping-beam intensities we 
studied. The barrier continues to operate asymmetri- 
cally, although the efficiency is strongly affected. 

The transmitting-side data in Fig. 10 show the ex- 
pected weak optimum in repumping-barrier-beam inten- 
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FIG. 7: (Color online) Populations on the left- and right-hand sides of the one-way barrier after 100 ms as functions of barrier- 
beam separation. From left to right, top to bottom, atoms were initially loaded on the left in the F=l ground state, on the left 
in the F=2 ground state, on the right in the F—1 ground state, and on the right in the F=2 ground state. Error bars indicate 
statistical error from at least 20 repetitions. 



sity. When the repumping barrier beam is too weak to 
pump atoms efficiently to the reflecting state, the atoms 
continue to transmit through the barrier, resulting in 
very little asymmetry, save for whatever side the atoms 
happen to be on after 100 ms. However, even at the 
weakest intensity used in Fig. 10, the repumping barrier 
beam still pumped atoms to the reflecting state and the 
barrier acted asymmetrically, although less efhciently. As 
the intensity of the repumping barrier beam is increased, 
the tails of its beam profile, which extend to the left 
(transmitting side) of the main barrier beam, will even- 
tually have enough intensity to optically pump atoms on 
the left-hand side of the barrier into the reflecting state 
before they reach the main barrier beam. In this limit 
the barrier becomes reflective on both sides, such that all 
atoms will remain to the left of the barrier. Though we do 
not reach this limit in our experiment, these two limits to- 
gether indicate the existence of an optimum repumping- 
barrier-beam intensity that we do observe in Fig. 10. 

For the reflecting-side data in Fig. 10, we see the same 
decrease in effectiveness as the repumping barrier beam 
becomes too weak to pump atoms to the reflecting state. 
However, we do not see the same decrease in barrier ef- 



flciency as the repumping-barrier-beam intensity is in- 
creased, since a stronger repumping barrier beam just 
causes the barrier to reflect better. These behaviors 
in the high- and low-intensity limits are evident in the 
reflecting-side data shown in Fig. 10, where the barrier 
ceases to effectively block atoms as the intensity decreases 
but its performance improves and then saturates (pre- 
sumably at the point where all atoms reflect) as the in- 
tensity increases. At i = 0, the atoms in the dipole trap 
have just enough spread that the tails of the distribution 
cross the barrier, as can be seen at t = in Fig. 11. This 
is a likely explanation as to why the population on the 
left-hand side of the barrier is not zero in Fig. 10, even 
for the largest repumping-barrier-beam intensities. 

We measured population losses for atoms released 
on the right-hand side of the barrier for two different 
main-barrier-beam detunings (shown in Fig. 12). Fig- 
ure 13 shows the results. The near-detuned and far- 
detuned data were collected with the main barrier beam 
tuned 0.75(5) GHz (using a Fabry-Perot cavity) and 
3.97(7) GHz (using the ^^Rb F=2 transition) 
blue of the ^^Rb MOT trapping transition, respectively. 
The near-detuned value was picked to be substantially 
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hand sides ol the one-way barrier after 100 ms as functions 
of initial loading position. The initial loading position is the 
position of the zero of the MOT magnetic fields relative to 
the focus of the dipole trap, in millimeters. Our canonical 
starting locations were about 0.65(8) mm away from the focus 
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regions. Vertical error bars indicate statistical error from 38 
repetitions, and horizontal error bars indicate uncertainty in 
measuring the MOT center. 
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FIG. 9: (Color online) Some sample time-series that formed 
the data in Fig. 8, showing populations on the right- and 
left-hand sides of the barrier as functions of time. Shown 
are curves for initial loading positions of (a) 0.89(8) mm, 
(b) 0.43(8) mm, (c) 0.14(8) mm, (d) -0.51(8) mm, and (e) 
—0.90(8) mm from the barrier center. Error bars indicate 
statistical error from 38 repetitions. 



closer to the F=2 F' resonances than our canoni- 
cal value, but not so close as to increase scattering of 
the barrier-beam light enough to prevent barrier oper- 
ation. The far-detuned value was conveniently close to 
halfway between the resonances for F—1 (transmitting) 
and F=2 (reflecting) atoms, and is actually a little closer 
to the F=l resonances. This allowed us to compare 
a more symmetric detuning to our canonical case and 
test how damaging state-changing scattering events are 
to the barrier operation. For both of these data sets. 



the intensity of the main barrier beam was adjusted so 
that the height of the reflecting barrier was the same 
as for our canonical data (with the barrier tuned to the 
^^Rb MOT transition), which is also shown in Fig. 1.3. 
The near-detuned data are nearly indistinguishable from 
our canonical reflecting-side data. The far-detuned data 
show a much shorter lifetime, but no extra leakage to the 
reflecting side of the barrier. This indicates the following: 

1 . Heating from scattering barrier-beam light alone is 
not much of a problem. If it were, then the near- 
detuned data, which should have ~ 2 times as much 
scattering as our canonical data because it is closer 
to resonance, should show much more loss. 

2. State-changing from scattering barrier-beam light 
is an issue. We believe we can explain this as fol- 
lows. In both our canonical data and the near- 
detuned data, the barrier was much more likely 
to pump atoms to the transmitting state, which 
decreases the potential energy of the atoms, thus 
helping to cool them. The repumping barrier beam 
pumps to the reflecting state, but since that oc- 
curs away from the barrier beam, heating due to 
the potential-energy increase is small. In the far- 
detuned data, the main barrier beam can pump 
atoms to either state. Atoms in the transmitting 
state that are pumped to the reflecting state are 
likely to be ejected from the barrier at high speed. 
Atoms in the reflecting state that are pumped to 
the transmitting state are more likely to be mov- 
ing slowly, because they were being reflected, which 
makes them likely to scatter again. Thus, atoms 
leaving the beam after a state change are more 
likely to have been heated than cooled. This in- 
creased heating decreases the trap lifetime substan- 
tially. Appendix B discusses scattering in more 
depth. 

We also briefly tried having both the main barrier 
beam and the repumping barrier beam polarized parallel 
to the longitudinal axis of the dipole-trap beam. We saw 
no significant change, suggesting that the barrier opera- 
tion is insensitive to polarization. 



V. SIMULATIONS 

We simulated the barrier operation using a simple 
model where atoms were assumed to be noninteracting 
point particles with well-defined positions and momenta. 
These atoms move in conservative potentials formed by 
the dipole-trap beam and the main barrier beam. At 
every point in time the atoms were assumed to be in ei- 
ther the F=l (transmitting) or F=2 (reflecting) ground 
state. This is because coherences between these states os- 
cillate at much faster time scales (subnanosecond) than 
the center-of-mass motion (microseconds to milliseconds) 
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FIG. 10: (Color online) Populations on the left- and right-hand sides of the one-way barrier after 100 ms as functions of 
repumping-barrier-beam power. The left panel shows results where the atoms were initially on the left-hand (transmitting) 
side of the barrier, and the right panel shows results where the atoms were initially on the right-hand (reflecting) side. The 
vertical lines in the plots show our usual repumping-barrier-beam power of 0.36 /iW. The two horizontal bars show the ending 
populations from data taken with no barriers [the one shown in column (a) of Fig. 4 and a similar data set where atoms 
were initially on the other side of the barrier] . The no-barrier cases show unequal left- and right-hand populations due to the 
oscillations of the atoms. Error bars indicate statistical error from at least 38 repetitions. 
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FIG. 11: (Color online) Populations on the left- and right-hand sides of the one-way barrier for a few different repumping- 
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FIG. 12: (Color online) A schematic of the rubidium emis- 
sion spectrum around the D2 line, showing the detunings we 
tested for the main barrier laser. The horizontal axis mea- 
sures detuning from the *''Rb F=2 F'=3 MOT trapping 
transition. Both *^Rb and *^Rb resonances are shown. Our 
canonical detuning coincides with the *^Rb F=3 —* F'=3,4 
crossover transition, and the far-detuned value coincides with 
the ^'^Rb F^2 F'=2 transition. For the near-detuned 
value, we monitored the detuning with a Fabry-Perot cav- 
ity. Each detuning from the *^Rb MOT trapping transition 
is shown with a solid arrow. This approximates the detuning 
seen by atoms in the F=2 (reflecting) state. A dotted arrow 
shows the detuning from the *^Rb MOT repumping transi- 
tion, which approximates the detuning seen by atoms in the 
F—1 (transmitting) state. 
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FIG. 13: (Color online) Populations on the left- and right- 
hand sides of the one-way barrier, comparing barrier per- 
formances when the main barrier laser operates at differ- 
ent frequencies. For the canonical, near-detuned, and far- 
detuned data, the main-barrier detunings were 1.05(5) GHz, 
0.75(5) GHz, and 3.97(7) GHz blue of the *^Rb MOT transi- 
tion, respectively. See Fig. 12 for details. The near-detuned 
and canonical data are very similar, while the far-detuned 
data show much more atom loss. Error bars indicate statisti- 
cal error from at least 28 repetitions. 



of the atoms, and so can be time averaged, effectively al- 
lowing each atom to be in a classical mixture of the two 
ground states. 

The potential the atoms move in is a sum of potentials 
from the dipole-trap beam and the main barrier beam. 
We assume the strengths of the potentials to be propor- 
tional to the local beam intensities. We use far-detuned 
approximations for computing the potentials, and we also 
use the rotating-wave and two-level-atom approximations 
for the barrier beams. The repumping barrier beam 
pumps atoms out of the F—1 ground state very quickly, 
and is too weak and too far detuned to produce a strong 
potential for atoms in the F=2 ground state, so we ignore 
its potential. After each time step of the integration, we 
compute a local scattering rate for each atom given the 
beam intensities at the atom's position and the current 
atomic state. The dipole-trap beam is very far detuned 
and has a scattering rate on the order of 3 s~^, which we 
ignore — in our simulation only the barrier beams scatter 
off the atoms. We then randomly decide whether each 
atom scatters a photon in this time step, with a prob- 
ability equal to the scattering rate times the time step. 
A scattering event is modeled by a randomly directed 
(with a dipole-emission distribution) single-photon recoil 
kick applied to the momentum, and a randomly chosen 
atomic state change based on the probabilities computed 
in Appendix A. 

We found that the simulations were somewhat sensitive 
to initial conditions (compare Fig. 14 with Fig. 15). We 
achieved very good agreement between simulation and 
experiment with the following initial conditions: Atoms 
were initially placed in an oblong Gaussian ellipsoid, 
placed to match the initial conditions of the atoms in 
the experiment we were modeling. The long axis lies 
along the dipole-trap axis, with a length comparable to 
that measured in our experiments (we used a standard 
deviation of 400 fiia). The two short-axis widths were 
chosen so that most of the atoms started out trapped 
in the dipole trap, with some small variations to help 
match the simulations to the data. The momenta for each 
atom were initialized to random values with a ~100 fiK 
Maxwell-Boltzmann distribution. We also simulated the 
nonzero loading time. For each atom we picked a uni- 
formly distributed random start time that was less than 
the loading time. Each atom is frozen until its chosen 
start time, at which point its position and momentum 
begin to change according to the rules described above. 
This emulates the actual loading procedure, where atoms 
enter the dipole-trap beam and become trapped through- 
out the loading time. Atoms just entering the trapping 
region do so with random momenta because atoms in a 
MOT scatter often. Once trapped in the dipole trap they 
become dark to the MOT beams and begin evolving ac- 
cording to the dipole-trap potential. Thus, atoms that 
have spent time in the trap have a preferred direction and 
position different from atoms that have just been loaded, 
as they have begun moving toward the dipole-trap focus. 
This gives a slight position-momentum correlation that 
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just-loaded atoms do not have. The intention of having a 
load time in the simulation was to model this correlation. 

We also noticed that atoms loaded differently on one 
side of the barrier than on the other. This is likely due 
to fringes in the MOT trapping beams. Thus, in our 
simulations, we used two slightly different sets of initial 
conditions. Atoms loaded on the transmitting side of the 
barrier had a temperature of 100 /iK, with a transverse 
standard deviation of 16 /xm. Atoms loaded on the re- 
flecting side of the barrier had a temperature of 80 fiK, 
with a transverse standard deviation of 8 /xm. In both 
cases, the atomic positions were centered about a point 
on the trap axis, 0.9 mm from the dipole-trap focus and 
the barrier beams. 

Figure 14 was derived from the same data as columns 
(b) and (d) of Fig. 4, but with simulation results added 
in for comparison. The dashed simulation curve shows 
the results of the simulations as described so far. As an 
example of how the simulations are sensitive to initial 
conditions. Fig. 15 shows the same experimental data 
as Fig. 14 compared to different simulations where the 
initial conditions had the same longitudinal spreads but 
zero transverse position and momentum spreads. These 
initial conditions have less kinetic energy than those used 
in Fig. 14, yet atoms pass through the barrier much more 
effectively. This is caused by an orbital effect in the 
dipole-trap beam, where atoms with angular momentum 
are slowed in their progression toward the trap focus (see 
Appendix C). As can be seen by comparing Fig. 14 with 
Fig. 15, the presence of angular momentum (Fig. 14) does 
effectively slow the atoms. In some of our simulations, 
this slowing accounts for some of the atoms remaining 
on the transmitting side of the barrier after 100 ms, be- 
cause they approach the barrier too slowly to even reach 
it on that time scale. This effect can also be seen in our 
experimental data (i.e.. Fig. 5). 

Figure 16 shows data taken from an experiment much 
like the one in Fig. 5, except with a slightly lower barrier- 
beam power and carried out over much longer time scales. 
The simulations as described (dashed line) completely 
miss the experiment — we needed to include light-assisted 
coUisional losses to best fit the data. In cold-atom traps 
there are several mechanisms where light affects colli- 
sions between atoms [27, 28, 29, 30, 31], which we collec- 
tively refer to as light-assisted collisions. In all of these 
an incident laser beam excites an atom which then de- 
cays. While the atom is excited, it either is affected by 
a different atom-atom potential or decays to a different 
state through an interaction with another atom, result- 
ing in a large change in kinetic energy that ejects the 
atom (and maybe another) from the trap. Thus, all of 
these loss mechanisms may be collectively modeled as a 
density-dependent loss term with a coefficient that is pro- 
portional to the probability that atoms may be excited, 
although we believe that one in particular (radiative es- 
cape [31]) is responsible for our losses. 

We modeled light-assisted collisions by estimating the 
density of atoms in the barrier region and randomly re- 



moving them with a loss rate proportional to that density 
(and the beam intensity). More precisely, if we let N be 
the number of atoms in a small volume V, then we ex- 
pect the loss equation for dN/dt to contain a term of the 
form cIN'^ /V A"^ , where / is the local laser intensity, and 
A is the detuning of the beam. Then //A^ is approxi- 
mately proportional to the scattering rate of the atoms 
and the c coefficient is independent of volume. This is 
the loss term we used to model light-assisted collisions. 
We picked c ~ 3 x 10^ cm^ mW~^ s~'^ to match the data 
shown in Fig. 16, and used that value in all other simu- 
lations. 

The difference between simulations with and without 
light-assisted collisions can be seen in Fig. 16, where the 
simulation without light-assisted collisions is a very poor 
fit for the data, but the simulation with light-assisted 
collisions can be made to fit the data very well. The dif- 
ferences can also be seen in Fig. 14. They are not as pro- 
nounced on the shorter times scales, but the simulations 
with light-assisted collisions still model the data better. 
This suggests that light-assisted collisions are the dom- 
inant loss mechanism for our one-way barrier. Density- 
independent effects provide a decent (exponential) fit to 
the data in Fig. 16, but we know of no mechanism that 
could explain such a large loss rate. We can rule out any 
effect not directly related to the barrier beams, such as 
collisions with background atoms, due to the much longer 
lifetime of atoms in the dipole trap (about 20 s) without 
the barrier beams. 

If we divide our c coefficient by the appropriate 
detuning (squared) for the experiment done by Kup- 
pens et al. [31], we find a coefficient of about 2 x 
10~ cm mW"^s-^ This should compare to the K 
value of 1.1(5) x 10"!° cm^ mW"^ s'^ quoted in Ref. [31]. 
Our density-dependent loss simulations were too sim- 
plistic and approximate to expect better than order-of- 
magnitude agreement. 



VI. CONNECTION TO MAXWELL'S DEMON 

In his 1871 Theory of Heat, Maxwell [23] contemplated 
a hypothetical creature who, being able to distinguish in- 
dividual molecules in a gas, could separate hot molecules 
from cold molecules, resulting in a heat flow from a cold 
temperature to a hot temperature. This creature became 
known as Maxwell's demon, and its apparent ability to 
violate the second law of thermodynamics remained an 
unsolved quandary for decades. Everyday experiences, 
and even the thermodynamical definition of temperature, 
dictate that heat flows from high to low temperatures in 
an irreversible fashion. This implies that separating hot 
and cold must decrease the entropy of a system. One can 
also see how such a demon could be used to perform use- 
ful work by creating a temperature differential and using 
an engine to extract energy from that. 

Many variations of Maxwell's demon have been pro- 
posed, including a simpler demon that allowed particles 
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FIG. 14: (Color online) Populations on the left- and right-hand sides of the one-way barrier, comparing data (points) to 
simulations (solid and dashed lines). The solid simulation curve has an extra light-assisted collisional loss mechanism. 
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FIG. 15: (Color online) Populations on the left- and right-hand sides of the one-way barrier, comparing data (points) to 
simulations (solid and dashed lines). These simulations and data are the same as in Fig. 14, except that the simulation initial 
conditions had no transverse spread in either position or momentum. These simulations do not match the data as well as the 
simulations in Fig. 14, illustrating how basic barrier results, such as the time atoms take to reach the barrier and whether atoms 
oscillate across the barrier, are strongly affected by the initial conditions. The solid simulation curve has an extra light-assisted 
collisional loss mechanism. 



into part of a larger chamber, but not back out again — a 
one-way barrier [24]. Compressing the physical volume 
of a gas without changing the temperature also decreases 
entropy, and again allows for work extraction from a heat 
reservoir without the need for heat flow into a colder 
reservoir. 

Today, Maxwell's demon is usually handled through 
the concept of memory [24, 25]. In the case of the one- 
way barrier, the demon must somehow have a record 
of the atomic positions. If the demon knows the ex- 
act atomic positions, then to him each atom occupies a 
very tiny volume in phase space. By allowing the atoms 
into one side of the chamber, he is not really compress- 
ing phase space, but rather rearranging it. Rearranging 
phase space is possible via normal Hamiltonian evolution 
and is not a violation of the second law of thermodynam- 
ics. 

If the demon does not already know the positions, the 



demon must somehow monitor the system and determine 
when to allow atoms through. A measurement must take 
place, and some device must take action on the result. 
Therefore, the measurement result must be somehow im- 
printed in the device. For the next measurement, either 
the result must be recorded elsewhere, or the first result 
must be erased. If the result is erased, that is an irre- 
versible process which involves dumping entropy into the 
environment. If the result is not erased, then there is a 
memory record which goes from a known initial state to 
a disordered state, containing all the disorder that was 
removed from the atoms. Thus, entropy for the combined 
device system does not decrease. 

As a quick example, assume we have N addressable 
atoms distributed between two equal volumes Vi and V2 ■ 
A demon is preparing to make sure all the atoms are in 
V2, but the location of each atom within the two volumes 
is completely unknown. Say the demon measures the 
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FIG. 16: (Color online) Populations on the left- and right- 
hand sides of the one-way barrier, comparing data (points) 
to simulations (solid and dashed lines) for a long experiment 
showing the barrier lifetime. The trap was initially loaded 
nearly symmetrically. The dashed curve shows a simulation 
without including light-assisted collisions, which models the 
data very poorly. The solid curve shows the same simula- 
tions including light-assisted collisions, which matches the 
data very closely. For each curve (data and both simulations), 
the curves from top to bottom show the total population, the 
population on the right side of the barrier, and the population 
on the left side of the barrier. Error bars indicate statistical 
error from 19 repetitions. 



position of each atom one by one, and if they are found 
in Vi, then they are moved to V2. Once it is known 
in which volume an atom resides, a process that swaps 
containers need not involve an entropy change, as the two 
volumes are of equal size. The demon now has a record 
of which volume each atom was found in. In the simplest 
case, this memory is a series of N bits each set to 1 or 
2. We will assume this memory to be constrained by the 
second law of thermodynamics. We will also assume that 
the initial state of the memory is known, since writing a 
bit over an unknown state is the same as erasing, which, 
as described below, requires dumping entropy into the 
environment. 

If the demon erases the bit after each measurement, 
possibly to reuse it, then after each measurement the de- 
mon must reduce the phase-space volume of that bit by 
a factor of 2. This is because the erasure process must 
take two states (representing which volume contained the 
particle) and map them into one (the erased state). At- 
tempting to avoid this by measuring the state of the bit to 
decide how to erase it just transfers the problem to some 
other memory that needs to be erased. By the second law 
of thermodynamics, this entropy decrease of fee In 2 must 
be accompanied by a matching increase of entropy else- 
where, so that the total entropy does not decrease. After 
measuring and placing all N atoms, the demon must have 
dumped Nk^ In 2 of entropy to the outside environment. 

If the demon uses N different bits and does not erase 
them, then after the cycle is complete the memory is in 
one of 2^ equally likely states, since each atom was in 



either Vi or V2 with equal probability. To an outside ob- 
server the atoms are now in a smaller volume, but the 
state of the demon's memory has gone from a known 
initial state (zero entropy) to one of 2^ possible states, 
which is the same entropy increase as if the demon had 
erased the bits and dumped that entropy into the envi- 
ronment. This also exactly matches the entropy decrease 
of the atoms being compressed from two identical vol- 
umes into one. The observer cannot know the initial 
atomic distribution without having made measurements 
of his own, either of the atomic positions or the demon's 
memory. In each case there is an entropy transfer. 

Here is a short computer analogy. Assume the atoms 
are in a one-dimensional trap, with the position repre- 
sented by a number stored in memory. The volume is 
proportional to the number of states the number can rep- 
resent, and the entropy is proportional to the logarithm 
of that, which is proportional to the number of digits, or 
bits, in that number. The first bit of the number could 
represent whether the atom was located in Vi or V2 ■ The 
demon above could be thought of as taking atoms from 
Vi and placing them into the exact same place within 
V2. In this analogy, that means forcing the first bit of 
the number to 1, effectively decreasing the size of the 
number by one bit. The number of total digits (entropy) 
did not decrease; rather, one was split off and transferred 
to the demon's memory, and no longer represents part of 
the atom's position. 

This demon-powered barrier is similar to our one-way 
barrier. We can ignore extraneous effects such as photon 
recoil (which results in heating), unwanted scattering, 
and beam jitter, and still account for entropy. All that 
is needed is the irreversible step and the corresponding 
measurement record (memory) associated with it. The 
"memory" for our one-way barrier is the repumping bar- 
rier beam, and it accounts for the entropy even if we ide- 
alize it. The ideal repumping barrier beam changes the 
state of atoms from F=l to F=2 and leaves F=2 atoms 
alone. That irreversible step can take place because the 
transitions available to the F=l ground state are of a dif- 
ferent frequency than those from those available to the 
F—2 ground state. When an atom is changed from the 
F—1 state to the F=2 state, a repumping-barrier-beam 
photon is absorbed and a different photon is emitted. Re- 
pump photons cannot change the atom from F=2 back 
to F=l, but the lower-frequency emitted photons can. 
Completely ignoring random spontaneous emission direc- 
tions, we look only at whether the frequency of each pho- 
ton is resonant with the F=l — F' transitions or not. 
Each photon that is not resonant means that the number 
of trapped atoms was increased by one. 

Let the trapping volume occupy a fraction r of the to- 
tal volume V. It can be shown that the optimum setup 
is to have the largest optical depth possible, and as many 
untrapped atoms in the trapping region as possible. The 
number of untrapped atoms would normally be propor- 
tional to the volume ratio r of the trapping region, but 
could be lower if many atoms had just been trapped. In 
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this limit each repumping-barrier-beam photon increases 
the number of trapped atoms with probabihty 



riVi 



rNi + N2 



(1) 



where A^i is the number of untrapped atoms, rNi is the 
number of untrapped atoms in the trapping region, and 
N2 is the number of trapped atoms. 

Say we trap some fraction / of the atoms, and it 
takes P photons to do so. We can use the entropy of a 
two-part ideal gas to write out the change in entropy of 
the atoms (two-part because the trapped and untrapped 
atoms must be distinguishable if the barrier is to tell 
them apart). The atoms change from N atoms in a vol- 
ume y to {1 — f) N untrapped atoms in a volume V and 
fN trapped atoms in a volume rV. The change in en- 
tropy for the atoms is 

A5a = k^N [/Inr - (1 - /) ln(l - /) - /In/] . (2) 

Likewise, the repumping photons change from P identical 
photons to fN photons in one state and P — fN photons 
in a different state, spread throughout the same volume. 
However, that change is more likely to occur toward the 
beginning, when there are fewer trapped atoms that can 
become untrapped [Vtmp is larger when N2 is smaller in 
Eq. (1)]. We can approximate the entropy increase for the 
repumping beam by summing up the entropy increase of 
the beam for each atom trapped. Each term in the sum 
is the expected number of photons it will take to trap 
the next atom times the entropy of a single photon in a 
mixed state, which reflects the uncertainty of whether or 
not the photon was emitted by a newly trapped atom. 
The change in entropy of the repumping beam can thus 
be written as 

fN-l 

^s, = - ^ ii^^" + (1 - In (1 - ^«)] 



n=Q 

rn 



rn + {N — ?i) ' 



(3) 



where Vn is the probability of trapping another atom 
after n atoms have been trapped, as given by Eq. (1). 

Equation (2) shows a decrease in entropy from the de- 
crease in volume, but an increase in entropy due to the 
multiple atom types. The increase disappears in the limit 
of all atoms being trapped (/ — > 1). The repumping bar- 
rier beam shows a net increase in entropy due to the 
change in states of the photons. This repumping-barrier- 
beam entropy increase can be shown (at least in a large- 
number approximation) to always be at least as large as 
the decrease in atom entropy, so the net entropy change 
is non-negative. Thus, the changes in photon frequency 
are sufficient to uphold the second law of thermodynam- 
ics. This analysis does not include momentum diffusion 
for either atoms or photons resulting from photon recoil, 
so entropy actually increases more than is shown here. 

We created a version of such a demon. Figure 5 shows 
data from an experiment where we had atoms spread 



throughout a volume, and then activated a one-way bar- 
rier in the center. The transfer was not completely 
efficient — not all atoms were trapped on the right-hand 
side of the barrier, but the spatial compression is quite 
good. However, the spatial compression is almost com- 
pletely countered by heating from spontaneous scatter- 
ing of the barrier-beam light. We found that the over- 
all phase-space volume occupied by atoms during this 
experiment decreased by 7(2)%. Better compression is 
likely achievable should we optimize the setup for phase- 
space compression; presently we optimize the barrier for 
transmission rather than compression. To achieve bet- 
ter phase-space compression, we would use a smaller, 
slowly moving barrier starting initially off to one side 
of the atom distribution, as described by Ruschhaupt 
et al. [15]. At first, this barrier would capture only the 
atoms with the highest energies, and would do so near 
their turning points so that they had very little kinetic 
energy (allowing for a weaker barrier and less scattering). 
We would then adiabatically sweep the barrier along the 
dipole-trap axis, through the focus of the dipole trap. At 
each point, the barrier would trap new atoms near their 
turning points, where they had very little kinetic energy. 
The atoms that were already trapped would not gain ki- 
netic energy. Thus when the barrier passed through the 
dipole-trap focus, all the atoms would be left in the center 
of the trap, with less kinetic energy. 



VII. CLOSING 

In summary, we have demonstrated an all-optical 
asymmetric potential barrier capable of increasing the 
overall phase-space density of a sample of neutral alkali 
atoms. Furthermore, we have found that the barrier is 
robust, with fairly substantial changes having little ef- 
fect, including mechanical changes, optical changes, and 
changes in atom-state preparation. We have also devel- 
oped some theory and simulations describing the barrier. 

The authors would like to thank Eryn Cook and Paul 
Martin for their helpful comments. This research was 
supported by the National Science Foundation, under 
Project No. PHY-0547926. 



APPENDIX A: STATE CHANGING 
PROBABILITIES 

Here we review the computation of approximate prob- 
abilities for linearly polarized light to change the state of 
an atom, for use in our simulations. In the dipole approx- 
imation, an atom must be excited by the incident light 
field and then spontaneously decay to a different ground 
state in order for the state to change. We also assume 
that the incident light field is a single monochromatic 
beam. Therefore, we wish to compute the amplitude for 
an atom changing from one state (i^o,mo) to another 
(Fijmi), via an excited state (f ,m'), and then sum 
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over the excited states and square that sum to obtain a 
probability for changing states. We make the approxi- 
mation that the D2 multiplet is dominant, and further 
that the splitting between the lines is small compared to 
the barrier-beam detuning. With these approximations 
we find the probability to be 



(Al) 



E 



{Fumi\d,, \F',m') {F',m'\dq \Fo,mo] 



where q represents the polarization of the incident light, 
q' represents the polarization of the spontaneously emit- 
ted light (—1, 0, and +1 for the z component of an- 
gular momentum), and dq and dqi are the appropriate 
spherical-basis dipole operators. The polarization of the 
incident laser field determines q, but we sum over all 
q' and all intermediate states {F' , m!) within the subset 
of excited states that are close to the laser frequency. 
The q' sum is done after squaring the amplitude be- 
cause different emitted polarizations are distinguishable. 
The (F',m') sum is performed with the amplitudes be- 
cause the near-degeneracy assumption means the dif- 
ferent emissions are indistinguishable. Since this does 
not represent all possible excited states (or any ground 
states), the sum of \F'm') {F'm'\ is not the identity. 

Next, we average over the initial magnetic sublevel toq 
and the emission polarization q' because we do not mea- 
sure or simulate the magnetic sublevels of the atoms. 
This is also an approximation, and is partially justified by 
the fact that the barrier operation was not significantly 
affected by the incident light polarization. 

Finally, we can conserve angular momentum by requir- 
ing m' + q' = mi and toq + 9 = ■ We can then evalu- 
ate Eq. (Al) and perform the averages over the magnetic 
sublevels of *^Rb. Given that we use linearly polarized 
light {q = 0) in our experiment, we find 



-^bari'icr,i^— 1- 
-Pbarricr,i^— 2- 



*F=2 = 5/18 = 0.27 
.F=i = 1/6 = 0.16. 



(A2) 



These are the values we use for the probability of a 
barrier-beam scattering event changing the state of an 
atom. The repumping barrier light is tuned close to the 
F=l F'—2 transition, so we require F'=2 in the F' 
sum, with the following result: 



^ropump,F=l^F=2 — 1/2 — 0.5. 



(A3) 



We found that the simulations were not particularly sen- 
sitive to these probabilities, although the barrier per- 
formed better in simulations with these ratios than in 
simulations where all ratios were 50%. As a side note, 
we also computed the probabilities that took into account 
the different dctunings of the excited states by weighting 
the excitation amplitudes in Eq. (Al) by factors inversely 
proportional to the detunings. The probabilities then be- 
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FIG. 17: (Color online) Plot of Eq. (B7), showing the ex- 
pected number of scattering events for a two-level atom pass- 
ing through a Gaussian barrier. The barrier height k is de- 
fined in Eq. (B4), and the number of scattering events is given 
as A''' :— NwomvoT/h |A|, where A'' is the actual number of 
scattering events. 
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These changes were not enough to significantly alter the 
simulations. 



APPENDIX B: 



ANALYTIC SCATTERING 
SOLUTION 



We can compute how many photons, on average, a 
two-level atom scatters while crossing a Gaussian beam if 
we ignore momentum kicks due to spontaneous emission. 
The barrier width is small compared to the Rayleigh 
length of the dipole trap, so we ignore any changes in 
the underlying potential and treat an atom in free space 
crossing a Gaussian potential with a scattering rate pro- 
portional to that potential. For a two-level atom in the 
far-detuned, rotating- wave approximation, the scattering 
rate is i? = |ry/?iA|, where F is the transition linewidth, 
V is the ac Stark shift (optical potential), and A is the 
(angular) detuning. The average number of scattering 
events can be computed by integrating the local scatter- 
ing rate for all time, which gives 



R[x{t)] dt, 



(Bl) 



where x{t) is the location of the atom at time t. For a 
Gaussian potential, we know the atom will approach the 
barrier and then either refiect or transmit. The expected 
number of scattering events up until the refiection (or 
passage through the barrier center) is the same as after 
the event, so we can restrict the integral to a period where 
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FIG. 18: (Color online) Plot of Eq. (B7), showing the ex- 
pected number of scattering events for a quasi-two-level (ig- 
noring state changes) *^Rb atom passing through a Gaussian 
barrier as a function of beam detuning from the _F=2 F' 
transitions. The curves that increase for larger detunings are 
for atoms in the F=\ transmitting state, and the curves that 
decrease for larger detunings are for atoms in the F=2 re- 
flecting state. The effective detuning from the F=\ ground 
state for a beam tuned to the F=2 — > F' transitions is about 
6.53 GHz. Each F=\ and F=l pair of curves represents a 
different barrier height. The dotted curves show the expected 
numbers of scattering events for a barrier height of = 1.25. 
This means that the beam intensity for each detuning is ad- 
justed so that the potential barrier seen by atoms in the F=l 
ground state is 1.25 times the kinetic energy of the atoms. 
The short-dashed curves are for k — 1.5, the medium-dashed 
curves are for k = 2, and the long-dashed curves are for = 4. 
Our canonical data had k = 2.3 with an effective two-level de- 
tuning of 1.12 GHz. 



the velocity does not change sign. This allows the change 
of variables to x, using dt = dx/v{x), where v{x) is the 
velocity of the atom when it is at position x. After this 
change of variables, we find 



iV = 2 



R{x) 
v{x) 



dx. 



(B2) 



where a;o = in the case where the atom passes through 
the center of the barrier, and the turnaround point in the 
case where it does not. 

If we let be the 1 /e^ intensity radius of the Gauss- 
ian beam, then we can use conservation of energy to com- 
pute 



j{x) = vq\/i 



(B3) 
(B4) 



where vq is the initial velocity of the atom [vq = v(x 
oo)] , Vq is the height of the barrier, and k is defined as the 
ratio of the barrier height to the incident kinetic energy 
of the atom. We can compute xq by solving v{x) = 0: 



Xo 




k < 1 
k>l 



(B5) 



A negative value for k means the barrier is actually an 
attractive well, and A: < 1 means the barrier is not high 
enough to stop the atom. For these cases, Xq = 0. Sub- 
stituting the scattering rate R(x) and the velocity v{x) 
from Eq. (B3) into Eq. (B2) gives 



N 



mV()T \ k\ 



dx, 



(B6) 



where Xq is given by Eq. (B5). 

We now change the integration variable in Eq. (B6) to 
u — ex-p{—2x^/wQ), with the result 



N = 



womvoT\k\ fl'^'^^i^Ml 
2V2h\A\ 



du 



(B7) 



The k limits are as expected 
1 



iV — > as A; ^ 0. The number of scattering events, 
N , converges linearly to zero as the barrier height 
k (or well depth, if k is negative) decreases. When 
there is no barrier, there is no scattering. 

2.iV^ooasfc^l. This is the usual logarithmic 
divergence (it diverges as — In |1 — A: |) of the time it 
takes an atom to roll off a hilltop starting with no 
initial velocity. In this limit, the number of scatter- 
ing events, TV, diverges because the barrier height 
k is approaching the limit where the atom comes 
to rest at the top of the barrier (and stays there for 
an infinite amount of time). 

3.iV— s-ooasfc^— oo. The number of scatter- 
ing events, N, diverges as y^jfcl as the potential 
well depth |fc| becomes infinite, meaning that for a 
very deep attractive well the increase in speed with 
which the atom passes through the barrier is not 
enough to counter the increase in scattering rate. 
Interestingly, if we keep the well depth constant and 
let fc — > — oo by decreasing the initial kinetic en- 
ergy (let Vq oc l/y^|fc|), N converges to a constant. 
This is because in this limit, the initial speed of the 
atom is negligible compared to the speed increase 
as it crosses the attractive well, so the initial speed, 
and thus fc, does not matter. 

4. TV^Oas/c^oo. The number of scattering events 
N converges as 1/Vln k as the barrier height k be- 
comes infinite, and means that if the barrier is much 
higher than it needs to be to refiect an atom, the 
number of scattering events can be decreased to 
zero. 

Equation (B7) shows why changing how far back we 
initially start the atoms is not quite the same as changing 
the barrier height. Changing the barrier height affects 
only fc, whereas changing the initial location also affects 

Vq. 

Figure 17 shows a plot of Eq. (B7). Note the weak 
dependence on barrier height for a reflecting barrier of 
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at least twice the required height (fc > 2). Also note 
the relatively strong dependence of scattering events on 
well depth for attractive wells (fc < 0). This hints at an- 
other reason why we want the main barrier beam tuned 
closer to the F—2 ground-state transitions than the F—1 
ground-state transitions. The main barrier beam will 
not scatter much on reflection if it presents a reason- 
ably high barrier, even if the detuning is not very large. 
However, the transmitting potential well that the main 
barrier beam presents to the F=l ground-state atoms 
quite possibly will scatter many photons, because the 
scattering rate increases much faster with |fc| for an at- 
tractive well. This can be seen in Fig. 18, where we plot, 
for a few different barrier heights, the expected num- 
ber of scattering events for each state of the atom as a 
function of beam detuning. The main-barrier-beam in- 
tensity was constrained to maintain a constant (indepen- 
dent of detuning) barrier height as seen by F=2 atoms. 
Since atoms only need to transmit once but must re- 
flect many times, one might think from Fig. 18 that the 
optimum barrier frequency is around 4 GHz from the 
F=2 F' transitions. However, as discussed earlier, 
scattering events for transmitting atoms are much worse 
than for reflecting atoms because they may change the 
state of the atom to the reflecting state, causing much 
more heating and loss than the reverse situation. We 
found that a detuning that reduced scattering events on 
transmission much more than on reflection was optimum. 
For the experiments we performed (with the exception of 
experiments where we changed how far back the atoms 
started), k was about 2.3 for atoms in the reflecting state. 
For atoms in the transmitting state (with the exception 
of experiments where either the frequency of the barrier 
or the initial position was varied), k was about —0.47. 
With our effective two- level detuning value of 1.12 GHz 
for the reflecting barrier, this yields ^ 8 scattering events 
on reflection, and ^0.7 scattering events on transmission 
(with an effective two-level detuning of —5.41 GHz). 



APPENDIX C: ANGULAR MOMENTUM IN 
SINGLE-FOCUSED-BEAM DIPOLE TRAPS 

In this appendix we briefly describe how angular mo- 
mentum affects atomic motion in a dipole trap formed by 
a single focused Gaussian beam. The basic effect is that 
an atom with angular momentum about the trap axis 
feels an effective force that repels it from the trap focus. 
A similar effect is experienced by a ball rolling inside 
a funnel. The tighter confinement toward the center of 
the funnel (or dipole-trap focus) forces a higher angular 
velocity in order to conserve angular momentum, which 
results in an outward push. This effect is often dealt with 
in classical mechanics textbooks such as Ref. [32] when 
dealing with central potentials. The angular momentum 
of a particle in a central potential leads to a centrifugal 
potential-energy term that repels the particle from the 
center of the potential. 
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FIG. 19: (Color online) Plot of cylindrical coordinates for two 
orbits near the focus of a dipole trap using the same dipole- 
trap parameters as in Sec. V, except without gravity. Longi- 
tudinal distances («) are in millimeters, and radial distances 
(r) are in ten-micron units. The solid curves show a simu- 
lation with 90% of the critical angular momentum required 
for the dipole-trap focus to become repelling. This trajectory 
shows an oscillation about the focus that is slowed due to the 
angular momentum. The dashed curves show the same initial 
conditions, but with a slightly larger r value and azimuthal 
velocity, pushing the angular momentum to 110% of the criti- 
cal value, and show how such a particle can appear to bounce 
off an attractive potential. 
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FIG. 20: (Color online) Plots of cylindrical coordinates for 
three orbits in a dipole trap, using the same dipole-trap pa- 
rameters as in Sec. V, except without gravity. Longitudinal 
distances (z) are in millimeters, and radial distances (r) are 
in ten-micron units. The solid curve shows a slight variation 
of the predicted saddle-point orbit, with the correct angular 
momentum but with q slightly below unity. This curve shows 
the particle oscillating about the saddle-point orbit, but drift- 
ing along the q=\ manifold. The radial restoring force for this 
orbit in our trap is much stiffer than the longitudinal restoring 
force, so r oscillates rapidly (resulting in a noncircular orbit). 
The dashed curve shows the predicted saddle-point orbit, but 
with 101% of the saddle-point angular momentum; this orbit 
diverges. The dotted curve shows the predicted saddle-point 
orbit, but with 99% of the saddle-point angular momentum; 
this trajectory falls away from the orbit and instead oscillates 
about the focus. 
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For this appendix, we will assume we have only the 
focused Gaussian dipole-trap beam, creating a conserva- 
tive potential proportional to the local intensity. The 
potential can thus be written as 



exp 



2r' 



1 + 



1 + 



(CI) 



where r is the radial coordinate measured from the trap 
axis, scaled by the l/e^ intensity radius at the focus, and 
z is the longitudinal coordinate measured from the trap 
focus, scaled by the Rayleigh length of the beam. We 
are discussing attractive traps, so Vq is assumed to be 
positive. The quantity in the exponential in Eq. (CI) will 
occur several times in this appendix, so we will abbreviate 
it as q^: 



2r' 



q ■= 



1 + Z2 



(C2) 



With this abbreviation, the equations of motion for the 
cylindrical coordinates become 



d?z 



2K) /, 2- 
2^ (1 ^ 9 



exp (-g2) 

(1 + ^2)2 



(C3) 



[j^~mwlVoq^exp{-q^)] . (C4) 



Here J is the conserved angular momentum about the 
trap axis, and zq and wq are the Rayleigh length and l/e^ 
intensity radius, respectively. Not surprisingly, the focus 
(z=0) is a fixed point of the longitudinal motion. What is 
more surprising is the {1 — q^) factor. The presence of this 
factor means that circular orbits with large angular mo- 
menta can have significantly smaller accelerations along 
the trap axis than for one-dimensional motion, even at- 
taining a different sign. We observed experimentally that 
our trap period was longer than what we would expect 
from a purely one-dimensional treatment. In simulations 
we could see that this was partially due to atoms in high- 
angular-momcntum orbits taking a longer time to reach 



the trap focus. We could also observe the repulsion from 
the trap center in simulations. This occurs when the en- 
ergy required for a given angular momentum at the focus 
is larger than the trap depth. Figure 19 shows simula- 
tions of some trajectories with these effects. We took the 
same simulations used in Sec. V, but with single atoms 
and initial conditions appropriate for orbits (and no grav- 
ity, in order to preserve cylindrical symmetry) near the 
critical point, illustrating how the focus of the dipole trap 
changes from attracting to repelling. 

This effect helps explain the difference between Figs. 14 
and 15. Atoms with some angular momentum have q > 0, 
reducing the longitudinal acceleration toward the trap fo- 
cus in Eq. (C3). This is why the simulation in Fig. 15 
shows atoms crossing the trap center sooner than the sim- 
ulation in Fig. 14. Furthermore, because atoms without 
angular momentum experience more longitudinal accel- 
eration, they are traveling faster along the trap axis when 
they cross the barrier than their counterparts with angu- 
lar momentum. We believe this is the main reason why 
the atoms traverse back and forth across the barrier in 
Fig. 15 more easily than in Fig. 14 — the barrier is ef- 
fectively higher for atoms that are moving more slowly 
along the trap axis. 

Linear stability analysis of Eqs. (C3) and (C4) sug- 
gests that for q^ ~ 1, there exist "saddle-point orbits" 
at some point that is not the trap focus. These orbits 
require q^ — 1 and = mwQVo exp(— 1), and are sta- 
ble to perturbations in q but not to perturbations in J. 
They have no linear restoring force along the q^=l man- 
ifold. We were able to produce simulations that show 
these strange orbits, but they were so sensitive that we 
had to remove the effects of gravity from our simulations 
in order for them to persist. Figure 20 shows one of these 
orbits and demonstrates how, even with the right condi- 
tions, the orbit drifts along the (?^=1 curve. This drift is 
a result of nonlinear terms in the equations, which domi- 
nate because the linear term is zero. Figure 20 also shows 
how quickly orbits with the wrong angular momentum 
diverge. 
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